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This is the first of a series of three articles that treats fracture localization as a critical phenomenon. 
This first article establishes a statistical mechanics based on ensemble averages when fiuctuations 
through time play no role in defining the ensemble. Ensembles are obtained by dividing a huge rock 
sample into many mesoscopic volumes. Because rocks are a disordered collection of grains in cohesive 
contact, we expect that once shear strain is applied and cracks begin to arrive in the system, the 
mesoscopic volumes will have a wide distribution of different crack states. These mesoscopic volumes 
are the members of our ensembles. We determine the probability of observing a mesoscopic volume 
to be in a given crack state by maximizing Shannon's measure of the emergent crack disorder subject 
to constraints coming from the energy-balance of brittle fracture. The laws of thermodynamics, the 
partition function, and the quantification of temperature are obtained for such cracking systems. 
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I. INTRODUCTION 

When rocks and other disordered-solid materials are 
in compression and then have an additional deviatoric 
strain applied to them, small stable cracks irreversibly 
appear at random throughout the material. Each time 
the deviatoric strain is increased, more cracks appear. In 
the softening regime following peak stress, a sample will 
unstably fail along a plane localized at an angle relative 
to the principal-stress direction. We have accumulated 
evidence suggesting that such localization is a continuous 
phase transition. 

This is the first of three articles that develops a sta- 
tistical mechanics that allows the possible phase transi- 
tions in a cracking solid to be investigated. Many stud- 
ies have assumed that fracture is a thermally-activated 
process and have used a statistical mechanics based on 
thermal fluctuations [1-5]. However, our interest here is 
with "brittle fracture" in which cracks appear irreversibly 
and in which thermal fluctuations play no role. For this 
problem, the statistics of the fracture process is entirely 
due to the initial quenched disorder in the system. 

A considerable literature has developed for so-called 
"breakdown" phenomena in systems having quenched dis- 
order and zero temperature [6-23]. In particular, the 
burned-fuse [6-8], spring-network [9-11], and flber-bundle 
[12-17] analog models for fracture have all been shown to 
yield various types of scaling laws prior to the point of 
breakdown [18-23] . Our work is different in that we di- 
rectly treat the fracture problem (not an analog model of 
it) assuming that all of the statistics is due to quenched 
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disorder. We obtain the probability of emergent damage 
states by maximizing Shannon's entropy subject to ap- 
propriate constraints. This approach has recently been 
proven exact in the special case of flber bundles [|3 . 



The principal conclusion of our present theory is that 
at a critical-strain point, there is a continuous phase tran- 
sition from states where cracks are uniformly distributed 
to states where coherently oriented cracks are grouped 
into conjugate bands. Several facts justify classifying 
such band formation as a critical phenomenon. 



First, the localization of the cracks into bands spon- 
taneously breaks both the rotational and translational 
symmetry of the material even though our model Hamil- 
tonian preserves these same symmetries. The entropy of 
the material remains continuous and the ensemble of the 
most probable states becomes degenerate at the localiza- 
tion transition; i.e., prior to localization, the most proba- 
ble state is the intact state, while right at the transition, 
certain banded states acquire the same probability as the 
intact state. Further, an autocorrelation length associ- 
ated with the aspect ratio of the emergent crack bands 
diverges in the approach to the critical point. Unfortu- 
nately, quantitative laboratory measurements of how the 
bands of cracks coalesce and evolve in size and shape 
prior to the flnal localization point do not presently ex- 
ist. We speculate in the third article of this series on how 
such measurements might be performed. 



Our explanation of localization based on the physics of 
interacting cracks is distinct from the bifurcation analysis 
of Rudnicki and Rice [|2^ in which localization is a conse- 
quence of a proposed phenomenological elasto-plasticity 
law. Our work provides a method for obtaining such a 
plasticity law from the underlying physics. 



2 



II. THE PROBABILISTIC NATURE OF THE 
FRACTURE PROBLEM 

Rocks are a disordered collection of grains in cohe- 
sive contact. The grains have varying shapes and sizes 
with typical grain sizes in the range of 10-100 jim but 
sometimes considerably larger. The contacts between 
the grains are generally weaker than the grains them- 
selves and have strengths and geometries that vary from 
one contact to the next. When deviatoric {i.e., shear) 
strain is applied to a rock, grain contacts begin to break. 
In what follows, a broken grain contact will be called a 
"crack". Such a break is a stress-activated irreversible 
process. Once a grain contact is broken, there is no 
significant healing that occurs. Cracks are not arriving 
and disappearing due to thermal fluctuations. This fact 
makes our deflnition of statistical ensembles quite differ- 
ent from that in the usual application of statistical me- 
chanics to molecular systems as we now go on to discuss. 



A. Creating a statistical ensemble 

We imagine dividing a huge (formally infinite) system 
into mesoscopic volumes that will be called "mesovol- 
umes". Because the materials of interest here have a wide 
range of grain-scale disorder, many different crack states 
will emerge in the various mesovolumes once energy has 
been put into the system and cracking begins. These 
various mesovolumes and the crack states they contain 
comprise the ensembles in our theory. 

In order to be specific with our ideas, we now intro- 
duce a simple model of the initial disorder and emergent 
crack states. The purpose of this special model in the 
present paper is to motivate how ensembles are formed; 
however, the model Hamiltonian developed in Paper II 
will be based upon it. 

In the model, each mesovolume is divided into N iden- 
tical cells where a cell has dimensions on the order of a 
grain size and where is a large number such as 10^-° or 
more with D the system's dimension. In each cell, only a 
single grain contact is allowed to break. The local order 
parameter (expHcitly defined in Paper II) characterizes 
both the orientation and length of such a broken grain 
contact. In the present paper, an order-parameter de- 
scription is not yet necessary. Prior to breaking, all cells 
are assumed to have the same elastic moduli. 

The quenched disorder is in how the grain-contact 
breaking energy S{x) is distributed in the cells x of a 
mesovolume. We assume that only a fraction of the nomi- 
nal grain-contact area is actually cemented together, and 
that the degree of cementation from one contact to the 
next is random. Thus, the breaking energies £{x) are 
random variables independently sampled from a distri- 
bution tt{£) having support on [0,rd^~^] where T is the 
surface-energy density of the mineral, d is the nominal 
linear dimension of a grain contact, and d^~^ is the grain- 
contact area in D dimensions. The quenched-disorder 



distribution Tr{£) can have any assumed form. 

We now define an infinite collection of distinct meso- 
volumes by allowing for every conceivable way that £{x) 
may be distributed in a mesovolume. Putting this collec- 
tion together forms the infinite rock mass whose proper- 
ties we are interested in determining. Each mesovolume 
so defined is a deterministic system and upon slowly ap- 
plying the same strain tensor e to all the mesovolumes, 
each will undergo a deterministic cracking scenario and 
end up in a well defined crack state. We denote each of 
the possible final crack states with an index j. A princi- 
pal goal of the present paper is to obtain the occupation 
probabilities pj of these various crack states which are 
simply the fraction of the mesovolumes in the system 
that are in the state j. 

We can understand how the various crack states 
emerge by appealing to a form of Griffith's crite- 
rion. A cell will break only if the change in the elastic 
energy due to the break is greater than or equal to the 
bond-breaking energy £{x). If Cq is the effective elastic- 
stiffness tensor of the entire mesovolume that holds after 
the break occurs and if Cfc is the stiffness-tensor that held 
before the break, Griffith's criterion can be stated 



'£:iCb-Ca):e/2>£ 



(1) 



where e is the strain tensor characterizing the entire 
mesovolume at the moment of the break and is the 
volume of a mesovolume. This particular statement is 
an approximation based on an assumed linear elasticity 
and absence of residual strain after unloa ding but a gen- 
eral statement will be derived in Section III B . Since the 



mesovolume with an extra crack is more compliant than 
without it, the weakest cells will begin to break even after 
the slightest of appHed strain. 

Yet an emergent crack state is not just a trivial conse- 
quence of the £{x) distribution in a mesovolume. Cracks 
aligned along bands concentrate stress allowing even 
large barriers £{x) to be overtaken along the band. In 
the present model, this means that placing cracks along 
bands produces a larger change in the elastic moduH of 
the mesovolume than placing cracks in more random po- 
sitions. Thus, at least above some applied strain level, we 
expect the banded states to emerge as the ones that are 
significantly present in a rock system. Non-banded states 
at large strain are much more special. They can come 
only from mesovolumes in which the weak cells making 
up the state are all surrounded by strong cells. 

A key idea here is that each mesovolume embedded in 
the system experiences the same global strain tensor and, 
as such, has a crack state statistically independent from 
the other mesovolumes. This is only valid so long as the 
emergent bands of organized cracks have a dimension ^ 
that is small relative to the size £ of the mesovolume. 
Screening effects due to destructive strain interactions 
between incoherently oriented cracks cause the far-field 
strain from a local crack structure to fall off with distance 
r even more rapidly than the {£,/r)^ fall off in an un- 
cracked material. But even in the thermodynamic limit 
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FIG. 1: Stress-strain data courtesy of David Lockner of the 
USGS Menlo Park. The slope measured upon loading a sam- 
ple is defined by D while that measured upon unloading 
and/or reloading the sample is defined by C. 



of infinite system sizes, the required statistical indepen- 
dence of the mesovolumes breaks down right at the crit- 
ical strain where divergent bands of cracks become im- 
portant. The conclusion is that although our ensemble- 
based statistics is valid in the approach to localization, it 
is incapable of describing the post-localization physics. 



B. Macroscopic observables 

In the laboratory experiments to which we apply our 
theory, a sample is immersed in a reservoir from which 
either uniform stress or strain conditions can be applied 
to the sample's exterior surface dfl. The macroscopic 
strain tensor e is defined in terms of the displacement u 
at points on dfl as 



1 



nu dS 



(2) 



where n is the outward normal to the sample's surface 
and is the volume of the sample in D dimensions. 
This definition of deformation thus corresponds to the 
volume average of the local deformation tensor Vu(a;) 
defined at interior points x of the sample. It will soon be 
shown to be conjugate to the macroscopic stress tensor r 
in the expression for the work carried out on the sample. 
If strain (rather than stress) is the control variable, the 
displacements at points x of the external surface dil are 
given by u = a; • e. 

As shown in Fig. ||, a typical compression experiment 
starts with the sample in a pure hydrostatic pressure 
state and then systematically increases the deformation 
in the axial direction, keeping the radial "confining" pres- 
sure Pc constant. Other ways of controlling the radial 
stress during the experiment are to keep a constant ratio 
between axial and radial stress, or to impose a constant 
radial deformation. So long as the confining pressure 
does not become so large as to induce a brittle-to-ductile 



transition these various experiments all result in the 
same type of localized structure at large axial strains. 
When axial strain is monotonically increased, cracks ar- 
rive at each strain increment and the deformation and 
stress changes are related as 



dT 



de 



■.de — D-.ds 



(3) 



where the fourth-order tensor D is called the tangent- 
stiffness tensor. This tensor defines the slopes between 
the various stress and strain components as the sample 
is being loaded and is an experimental observable. 

If at some point in the stress history the axial pressure 
is reduced, we follow a different deformation path as seen 
in the figure due to the fact that no new cracks are cre- 
ated. Such an unloading experiment defines the elastic 
(or secant) stiffness tensor C 



dr = C : de. 



(4) 



We model the unloading/reloading paths as being en- 
tirely reversible and in so doing neglect the small hys- 
teresis due to friction along the opened cracks. 

In order to distinguish loading paths (with crack cre- 
ation) from unloading paths (without crack creation) , all 
properties are explicitly taken to depend on two strain 
variables; namely, the maximum strain £m having been 
applied to a sample, and the current strain e that is differ- 
ent than the maximum only if the sample has been subse- 
quently unloaded. Note that even if e and e„i are written 
as tensors, they each correspond to only one scalar de- 
gree of freedom along the loading/unloading paths, since 
the radial components can always be expressed in terms 
of the axial components via the type of radial control 
employed (e.g., pc = const in a standard triaxial test). 

The stress tensor r in the theory corresponds to the 
volume average of the local stress tensor T{x) that sat- 
isfies V • T(a;) = at interior points x; i.e., t — 
/n T(a;) dV and is a function of the current and 
maximum strains t — T{e,em) as shown in Fig. 1. By 
averaging the elastostatic identity V • (T a;) = T over the 
mesovolume we further have that r = L^^ J^^ n-Ta; dS. 

The work density dU performed on the sample when 
there is an increment in strain de is in both cases of 
loading and unloading 



dU = 



1 

T:de 



n-T-dudS (5) 

(6) 

To obtain Eq. (^ from (||) , we have written the controlled 
displacements on a sample's surface as du — x ■ de where 
the strain increment de is uniform over dfl. Thus, dU 
corresponds to the volume average of the local work den- 
sity T{x):d\/u{x). 

The total energy U per unit sample volume that goes 
into the sample during the loading up to a maximum 
strain-tensor em is then 



C/(£,„) 



T{e',e'):de'. 



(7) 
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where £o is the strain associated with the initial isotropic 
stress. If after loading to e™, the sample is unloaded back 
to a current strain of e, we have the general expression 

U{£,e„,) ^U{e,n) + T{e',£,n)-d£'. (8) 

If the sample is unloaded back to the initial stress, cor- 
responding to a possibly non-zero residual strain e'-''^^, 
a last experimental observable is the energy Q{em) = 
U {e'^'^^, Em) (per unit sample volume) that went into crack 
creation and that is lost during the loading process . 



C. Ergodic hypothesis 

We have shown above that the experimentally- 
measurable variables of energy density U, deformation 
e, and applied stress t correspond to volume averages 
of each field throughout a system. Our ergodic hypothe- 
sis amounts to assuming that the systems we work with 
are sufficiently large that such volume averages can be 
replaced by ensemble averages 

U = ^PjEj, £ = ^pj£j, T^Y^PjTj. (9) 

j 3 3 

Here, Ej is the average work per unit mesovolume re- 
quired to take an initially uncracked mesovolume from 
zero strain to the strain tensor Ej. A similar defini- 
tion holds for Tj. In both the definition of Ej and 
Tj — dEj/d£j, the average is over the initial quenched- 
disorder distribution. 

So long as each mesovolume contains crack states that 
have no significant infiuence on the neighboring mesovol- 
umes (formally valid only in the thermodynamic limit), 
the sum over the collection of mesovolumes (ensemble 
averaging) is equivalent to a volume integral over the 
entire system. In practice, we will only ever consider 
ensembles that have by definition ej = e; however, we 
could equivantly immerse each mesovolume in a uniform 
stress-tensor reservoir and allow Sj to vary from state to 
state. 



III. THERMODYNAMICS OF CRACK 
POPULATIONS 

A. Fundamental postulate 

The fracture- mechanics problem of counting how many 
of the initial mesovolumes can be led to the same crack 
state appears to be hopelessly intractable. Fortunately, 
it also appears to be unnecessary for systems containing 
initial quenched disorder. Upon putting deviatoric strain 
energy into such a system, the emergent crack states j 
will on the one hand attempt to mirror this quenched 
disorder with weakest cells breaking first; however, due 



to the energetics of the crack interactions, many differ- 
ent types of initial mesovolumes may be led to the same 
crack state which results in non-uniform crack-state prob- 
abilities pj even if the quenched-disorder distribution is 
uniform. 

We state our fundamental postulate as follows: The 
probability pj of observing a mesovolume to be in crack 
state j can be determined by maximizing Shannon's Jj^ 
measure of disorder 

S^-J2P^lnp, (10) 

3 

subject to constraints involving the macroscopic observ- 
ables that derive from the energetics of the fracture me- 
chanics. That entropy is to be maximized can be ex- 
pected since the quenched disorder allows all states to be 
present in a sufficiently large system. In recent work |Q, 
we have demonstrated that this postulate yields exact re- 
sults for the special case of fiber bundles with global-load 
sharing. 

The constraints are what give the dimensionless func- 
tion S defined by Eq. ( [Tc| ) all the thermodynamic infor- 
mation about our cracking system and must explicitly 
involve the independent variables of S. Such indepen- 
dent variables are determined by establishing the first law 
of thermodynamics for a system cracking in compressive 
shear. 



B. The work of creating a crack state 

To obtain the first law, it is first necessary to define the 
detailed energy balance for each crack state and to under- 
stand how the work Ej required to create state j depends 
on both the actual strain e and on the maximum-achieved 
strain £m- 



1. Griffith's criterion and crack-state energy 

Consider a given mesovolume with a deterministic dis- 
tribution of breaking energies S(x) assigned to each cell 
X of the mesovolume. Starting from a state of isotropic 
strain £o, we slowly apply an additional axial deforma- 
tion and monitor how one crack after another enters the 
mesovolume until the final strain tensor e and final crack 
state j are arrived at. Lets say that this state j has a 
total of A/" cracks associated with it. 

Figure ^ details the history of how the stress (and, 
therefore, work) might evolve in the mesovolume as strain 
is applied and cracks arrive. Initially, the mesovolume 
will elastically deform according to the stiffness tensor 
Co (no cracks yet present) until the first crack arrives at 
the strain tensor ei with an associated drop in the meso- 
volume's stress. Lets say the bond-breaking energy of 
this first crack was £i . The mesovolume will now have a 
different overall stiffness tensor Ci and will elastically 
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FIG. 2: The heavy line is the actual path followed during the 
steady application of axial strain. Each vertical drop in stress 
corresponds to the arrival of a crack. 



deform with these new moduli until the second crack 
arrives and so on until all M cracks have entered and 
the mesovolume has attained its final stiffness tensor of 
Cj = Cj\f. The final tensor Cj depends on both the lo- 
cation and orientation of these N cracks in addition to 
their number. 

At some intermediate stage having n cracks, the stress 
tensor r„(£) is defined by integrating dr = C„(£') : de' 
from 6^°*^ to £, where e^^^ is the "residual" deformation 
observed upon unloading the sample back to zero stress 
as shown in the figure. We have 



T„(e) 



C„:de'. 



(11) 



The elastic energy density corresponding to this state at 
deformation e is similarly 



^{e'):de' 



(12) 



where E'lj^^ represents the residual elastic energy that re- 
mains in the system when the state with n cracks is un- 
loaded to zero applied stress. These residual (zero-stress) 
quantities are present whenever plastic deformation oc- 
curs within a grain-contact. After a sample elastically 
returns to zero apphed stress, such plastic deformation 
remains and, accordingly, there is an elastic stress field 
surrounding any crack that experienced plastic deforma- 
tion. The strain energy associated with such local resid- 
ual stress is what constitutes the residual energy E]^^. 

When the n**^ crack arrives in a strain-controlled ex- 
periment, there is no change in the strain e„ and thus no 
external work performed. However, there is a change 
in stiffness (and possibly residual strain) resulting in 
an associated stress drop At„ = T„_i(e„) — T„(e„), 
and a drop in the stored elastic energy density Ai5^' = 
E^_i{en) — E'^i^n)- Energy conservation requires the 
elastic energy reduction to exactly balance the work per- 
formed in opening the crack so that 



el 







(13) 



where £n is the bond-breaking work performed at the 
grain contact of the n**^ crack, Kn is the energy that went 
into acoustic emissions when the crack arrived and/or 
expended in any mode II frictional sliding or plastic de- 
formation at the grain contact {Kn is a positive "loss" 
term), and, as earher, £^ is the volume of a mesovolume. 
Because Kn is positive, we can rewrite Eq. ( |l^ ) as an 
inequality 



(14) 



which is a general statement of Griffith's criterion. Upon 
appeahng to linear elasticity (elastic stiffnesses indepen- 
dent of strain level) and putting the residual deformation 
to zero (no plasticity inside the cracks), we arrive at the 
convenient statement £^ Sn ■ (C„_i — C„) : £„/2 > f„ 
given earher. 

The work performed between the arrival of the n*-^ and 
the (n -I- 1)**^ crack is defined 



Wn 



Tnie'):e' - EfHsn+i) - E^\en). (15) 



Thus, the total work required to reach the final strain e 
is the sum (c./.. Figure ||) 



E^=j:wn 



(16) 



m— 



where by convention W^/■ is the work performed after the 
arrival of the last crack to get to the final deformation e. 
The superscript p on E^ is simply indicating that this is 
the work for one particular realization of the quenched 
disorder. Rewriting the sum by introducing Eqs. ( |l^ ) 
and (|l|), then gives 

AT 

E^ = EJ^{e)-Et{eo) + Y.AE';! 



n=l 



Efie) + T.^-'^-Eo\eo) (17) 



ri=l 



where Eq^so) is the small and physically-unimportant 
amount of energy that is stored in the initial isotropic 
strain field. Equation ( |l7| ) is the natural statement that 
the work performed in creating state j at strain e is the 
sum of the elastic energy density stored in the material 
in the final state plus the energy irreversibly expended 
during the opening of each crack. 

Both the loss term Kn and the residual energies E^"^^ 
(contained in i5|') are potentially a function of the point 
in strain history at which a grain contact actuahy breaks; 
e.g., most models one might propose for plastic defor- 
mation at a grain contact are dependent on the applied 
stress level. However, modeling such plastic processes 
seems uncertain at best. We thus assume that at least for 
those crack states significantly contributing to any phase 
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transition (states with lots of cracks), the stress-history 
dependence of Kn is, on average, negligible. Further, 
since the residual strain in brittle-fracture experiments is 
never more than a few percent of the peak-stress deforma- 
tion and since the essence of the localization process does 
not seem to lie in E^°^, we assume that Ej°^ <C J^n^^- 
With these approximations, the work density Ej depends 
only on the final state j, the final strain e (through E^^), 
and the breaking energies £„. 

The energy density Ej needed later in our probability 
law is obtained by further averaging over the quenched- 
disorder in the breaking energies £n to give 

E, = Ef{e) + 7,(£™)^ - Ef{e,). (18) 

Here, Nj = Af is the total number of cracks in state j and 
7j is the average energy required to break a single grain 
contact where the average is over all cells throughout 
all mesovolumes led to state j. This jj can be different 
for different final crack states. It will also be greater at 
greater values of the maximum strain Sm because, ac- 
cording to Griffith, the cells comprising j can break at 
higher energy levels when the strain is greater. The first 
term in Eq. ( [l^ ) corresponds to the purely reversible 
elastic energy and therefore depends only on the actual 
strain state £. 



where the fraction fj is bounded as < /j < 1. We next 
demonstrate that the variation of fj from one state to 
the next is so small as to be neglected altogether. 

A tighter lower bound for fj is obtained by considering 
crack states j having Afj non- interacting cracks. Since 
the cracks do not interact to concentrate stress, all of 
the Mj cells that broke had their breaking energies some- 
where in the range < £ < SE — i^Sm ■ SC : £m/2 where 
(5C is the change in the stiffness tensor due to the arrival 
of a single non-interacting crack and 6E is the associated 
change in the elastic energy. Since the breaking energies 
are independent random variables taken from the distri- 
bution 7r(f), we obtain 



7j = 



lo^ e7r(e) de 
la^ 7r(e) de 



(23) 



for non- interacting crack states j. 

We now appeal to a specific form for the probability 
distribution tt{£). Initially, our rocks are intact and it 
is expected that more grain contacts are entirely bonded 
{£ — Td^~^) than entirely unbonded {£ = 0). We thus 
assume a monotonic distribution £'' with fc > satisfying 



the normalization J^'^ n(e)de = 1 so that 



£ 



jk + l) 



c£\ 



(24) 



S. Specific expression for Ej 

To facilitate the development in Paper II and to be 
more specific, we now use Griffith's criterion to develop 
an expression for Ej that is based on linear elasticity. 
When the n*'^ crack arrives, the linear-elastic variant of 
the Griffith criterion gives that 

£n < £^£„:(C„_i -C„):£n/2 (19) 
< £^£™:(C„_i -C„):£m/2 (20) 

where as earlier e„ is the strain point on the load curve 
where the n^^ crack arrives while Sm is the final maximum 
strain level of the experiment. The second inequality fol- 
lows from the first since an extra crack always reduces 
the stiffness of a mesovolume. For any particular meso- 
volume in state j, the average energy required to break 
a contact 7^ thus satisfies 



^71 



^3 

n— 1 



(21) 



where the right-hand side comes from summing Eq. (20). 
Since this inequality is independent of the history, every 
mesovolume that is led to state j must satisfy it. We may 
thus write 7^ in the form 



Using this tt, the average energy required to break a con- 
tact in a non-interacting crack state is 



1 



k 



-5E 



■.SC: 



(25) 



where we have defined q — (fc-|-l)/(fc+2). All dependence 
on the underlying quenched-disorder distribution in our 
theory is confined to the constant q which for any fc > 
is in the range [0.5, 1]. 

Since for non- interacting states Co — Cj = AfjSC, a 
comparison of Eqs. (|25| ) and (|2^ ) shows that fj — q for 
all the non- interacting states. For the interacting states, 
the prefactor fj must be slightly greater because now 
stress concentration can allow stronger cells to break. It 



is thus concluded that for all states, the fj of Eq. ( |22| ) are 
bounded as q < fj < 1 which when compared to how Mj 
varies from state to state can be considered negligible. 
From here on, we simply take fj = q for all states. 

The essential physics for the average amount of work 
that goes into building up any given crack state j is thus 
captured by 



Ej{£,em) 



E^i 



E^{e)+E'^{em) 
2 ' 

— £„i : (Co — Cj): e,i 



(26) 
(27) 
(28) 



(22) 



where the superscripts R and / denote respectively the 
reversible and irreversible part of the energy. The intact 
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hydrostatic energy Eq^{£o) has been neglected since it 
does not involve cracks and, therefore, cannot influence 
the probability of the various crack states. 



C. The laws of our crack-based thermodynamics 

Using the ergodic hypothesis discussed earlier, the av- 
erage energy density in a disorded soHd can be written 
U = J2jPj^j- We are interested in how U changes when 
increments in £ and £,„ are applied to the system. 

In general, a small increment in U can be written 



dU Ej dpj + ^ dEj . 



(29) 



The flrst term involving the probability change is en- 
tirely due to crack creation. Some mesovolumes that 
were in less cracked states prior to the increment, are 
transformed to state j during the increment, while meso- 
volumes that were in state j, are transformed to other 
more cracked states. If in the increment, the number 
of mesovolumes arriving in state j is different than the 
number leaving, there is a change dpj in the occupational 
probability of that state. Such changes are the only way 
to change the disorder in the system, so that 



Ejdpj — TdS 



(30) 



The "first law" for the rock mass is then 



dU 



■.de + g: de.m + TdS 



(34) 



with the formal definitions 



dU 



and T 



'dS 



(35) 

The natural variables of the fundamental function U are 
(5, £, £,n). Equivalently if S" is treated as the fundamental 
function, then S = S{U,£,em) which means that the 
constraints placed on the maximization of S must involve 
U, £, and Em- 

The "second law" of this crack-based thermodynamics 
is that c?5 > (equal to zero only if dem = so that no 
cracks are created) while a "third law" may be proposed 
by simply defining T = when S = 0. The system 
will have zero emergent disorder before cracks begin to 
arrive and so our third law states that the temperature 
T starts at zero and then increases in magnitude as the 
number of cracks in the system increases from zero. The 
justification for this postulate comes a posteriori when 
it is found that in order to have zero probability for a 
mesovolume being in anything but the uncracked state 
{S = 0), we must have that T = 0. 



D. The probability distribution 



is the work involved in changing the system's disorder 
via crack production. The proportionality constant T is 
formally a temperature and will be treated in detail. 

Using the decomposition Ej{e, £„) = E^{e)+Ej (em), 
we can write the second term of Eq. ( p9| ) as 



To obtain the probability of observing a mesovolume 
to be in crack state j, we maximize S = —J^jPj^^Pj 
subject to the constraint that J2jPj = ^^'^ to the 
additional constraints that £,- — £, £„y = £ 



T.jPj 



and 

Ej = U . These constraints define our canonical 
ensemble. Other ensembles can be defined by consider- 



'y^PjdEj = "y^pjdEf + ^y^pjdE^j. (31) ing other constraints involving £, £m, and U; however 



The first part is due to purely elastic (reversible) changes 
in each mesovolume and may be further written 



^PjdEf ^T-.de 



(32) 



where r is the average stress tensor acting on the meso- 
volumes. This result can be verified by appealing either 
to Eq. ( p7|) or to the more general statement of Eq. (12). 

The second part PjdEj represents the average work 
performed in creating cracks in just the final strain incre- 
ment de„i. Some of the initial mesovolumes led to state 
j at maximum strain £m + dsm had all their cracks in 
place before the final strain increment, while others had 
cracks arrive in the final increment. We write 



'^PjdE' ^g:d£r. 



(33) 



where the tensor g has units of stress but is quite distinct 
from the stress tensor r. 



since all ensembles yield identical average properties in 
the thermodynamic limit, we elect to work only with the 
canonical ensemble due to its analytical convenience. 

This maximization problem is solved using Lagrange 
multipliers to obtain the Boltzmannian 



-E,/T 



Pj 



where Z^^e"^^/^, 



(36) 



and where the parameter T is exactly the partial deriva- 
tive dU/dS\e,£^ called "temperature". 



E. The free energy and its derivatives 

Any equihbrium physical property that depends on the 
distribution of cracks throughout the system can be ob- 
tained from the partition function Z given by Eq. (|3^). 

To do so, a thermodynamic potential F called the free- 
energy density is introduced that is related to Z by 



F(£,£™,T) = -rinZ(£,£„„T). 



(37) 
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This potential F is the Legendre transform with respect 
to S of the total-energy density U = U{e, £„, S) as can 
be seen from 

U-TS = J2p,E,+tJ2 Pi Inp, = -T In Zj^Pj^F 

3 j 3 

(38) 

where we used that Inp^ = —Ej/T — hiZ. 

When (e,em,r) are the independent variables, the 
first law can be obtained by taking the total derivative 
of Eq. @ 



dF 



_T— -InZdT 
Z 



dEj{e,er. 



T 



P3 



- In Z dT 



dT 



(F - U)— +Y.P= [dEf{e) + dEl {sr, 

3 

~SdT + T:d£ + g:d£,n 



(39) 



where we have used the definitions that Tj = dE^'{e)lde 

and = dE^{£,n)/d£„i- 

With (3 — 1/T, the various thermodynamic functions 
are related to the partial derivatives of In Z{£, £„, /?) as 



dliiZ 

^ dp 

1 ainZ 



P de 

1 dliiZ 



(3 d£n 



Y.^3P3=U 

3 

E '^3P3 = T 

3 

^93P3 =9- 



(40) 
(41) 
(42) 



These results, along with S — \nZ + PU, are used in 
Paper III. 



IV. THE TEMPERATURE 

The temperature is a well-defined essential part of our 
quenched-disorder statistics. Through the probability 
law pj = e~^^l^ jZ, the temperature quantifies the en- 
ergy scale that separates probable from improbable states 
and how this energy scale evolves with strain. No other 
meaning should be read into T. We now demonstrate 
how to exactly obtain T . 



A. Evolution of temperature with strain 

The only way energy enters the system is by perform- 
ing work on the external surface. Thus, the general rela- 
tion dU = T:d£ holds for either loading or unloading sit- 
uations. This previously unused fact provides a differen- 
tial equation for T = 1//3 that permits everything about 
our system to be exactly known once an order-parameter 



based model for Ej{£,£m) is determined and the func- 
tional sums defining Z(£, £„,/?) = e~^^j(^'^™) are 
performed. 

The temperature and entropy only evolve along load 
paths defined by e = £m and only such paths need be 
considered in what follows. Using dU = r : d£, the first 
law [Eq. (|^] can then be rewritten 



TdS + g:d£ = 0. 



(43) 



Since it always requires energy to break contacts, we have 
that g .d£ > and consequently TdS < 0. Furthermore, 
since the entropy (disorder) necessarily grows during the 
crack-creation process (at least initially), the tempera- 
ture of our system is negative (at least initially) . 

The load path of a standard triaxial experiment is 
when axial strain Ez monotonically increases while the 
radial confining stress Tx = Ty = —pc remains constant. 
Along this path, all properties evolve only as a function 
oi Ez- With Z(e,£m,/3) considered as known, the radial 
deformation components can be expressed in terms of the 
axial deformation by using the two equations 



Ppc^ 



dlnZ 



dlnZ 



de,j 



to obtain the two functions 

£x = fx{P,£z) and Sy 



fyiP,£z 



(44) 



that are valid only along the load path. 

We now write dU in two different ways. First, dU = 
T:d£ is evaluated along the load path to obtain 



dU = TzdSz - Pc{dfx +dfy). 

Second, we use the fact that U = U{P, £, £ 

'dU dU 
^ dsz de,n 
dU \ ,„ /dU dU 



dU 

dU^-dp- 



dsz 



dU 

dsx 



(45) 
to obtain 



dfy. (46) 



my 



Upon equating Eqs. ( p5| ) and ( |4q ) we obtain a first-order 
non-linear differential equation for P 



a{P,ez) 



dp 
dez 



■b{P,ez 







(47) 



Pc 



where a and b are given by 

dU 
''=dp 



b = -i 



Pc 

dU 
dez 



Pc 



dU 

+ 7^ + 

oex 


du ^ 


\ dfx 




1 dp 


dU 


du ^ 


dfy 


dey 




dp 


dU 






^ de 






dU 


du ^ 


1 dfx 


dex 




' dez 


dU 


du ^ 


dfy 


dey 




dez 



(48) 



(49) 
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We are to use = —(3~^dlnZ/dez and U = —dlnZ/dp 
in these expressions for a and b once the function 
Z{e,em,P) has been determined. Furthermore, all par- 
tial derivatives are to be evaluated along the load curve; 



fy{P,e^), and Sr, 



B. Initial conditions 

In order to integrate Eq. (|47|), initial conditions must 
be provided. The initial conditions of our so-called "third 
law" {i.e., the intact conditions that P = — oo when Sz = 
0) are not well-defined for (3. Thus, Eq. ( p7| ) must be 
integrated not from the intact state, but from a state 
that contains at least a few cracks so that /? ^ — oo. 

Accordingly, we define "one-crack" initial conditions by 
considering the point in strain history where on average 
throughout the ensemble of mesovolumes, there is one 
crack in each mesovolume. If there are N cells in a meso- 
volume, the probability of any given cell to be broken 
somewhere in the ensemble is then Pi = 1/N. This 
same probability can also be obtained from Griffith's 
criterion by integrating the quenched-disorder distribu- 
tion of Eq. (HI) to obtain Pi = [SEi/(rd^-'^)]''+^ where 
5Ei = i^ei :(5C:£i/2 is the elastic energy change due to 
a single isolated crack and where ei is the strain tensor at 
which on average there is a single crack in each mesovol- 
ume. Thus, we have ei:5C:£i = 2Td°~^ /{N^/'^'^+'^k") 
which can be used to obtain an expression for the initial 
axial strain Ezi at which on average there is one crack 
per mesovolume. 

To obtain the inverse temperature (3i corresponding 
to this initial strain, the exact probability of observing 
a particular type of crack state is determined and com- 
pared to our temperature-dependent Boltzmannian. The 
particular states we choose to analyze are, for simplicity, 
those having precisely one broken cell. 

The probability pj of a state consisting of one broken 
cell and iV — 1 unbroken cells can be written 



Pj = Pi{i-Piy'-'n^ [1-SPix)] 



(50) 



where Pi is again the probability of having a single bro- 
ken cell and (1 — Pi)^^^ is the probability of having TV— 1 
broken cells in the absence of other cracks. Thus, the 
product Tlx [1 ~ SP{x)] is the probability that no cells 
broke due to the strain perturbations caused by the pres- 
ence of a first broken cell where x represents distance 
from this first broken cell. We define 5E2{x) as the elas- 
tic energy change in a mesovolume when a second cell 
breaks solely in the perturbed strain field emanating from 
a first broken cell. This energy varies with the separation 
distance \x\ between the two cracks as |a;|~^. We have 



SP{x) 



SE2(a 



7r(e) de 



fe+1 



C2 



\x\D(k+l) 

_ (51) 
where Eq. (gj) was used for tt and where C2 depends on 
both the overall applied strain and the angle from the first 



crack's orientation to the second crack. Since SP is small 
compared to one (restricting to models where cracks are 
smaller than the cell size A, since the separation distance 
I a; I always exceeds it), we have 

Ux[l-SPix)]^l-l,J^^^^^^-^d-x (52) 

and since fc > 0, this spatial integral over the mesovolume 
can be neglected in the thermodynamic limit. 
The conclusion is that 



Pj=Pi(l-Pi)^-i=Po 



Pi 



1-Pi 



= Poe-'"(^^^) (53) 



where po — (1 — Pi)^ is the probability of the entirely 
intact state. This can be compared to our probability 
law where, from Eqs. (pG)-(ES), we have 



Pj = Po exp 



Pi £i :dC:£i 



(54) 



Thus, the inverse temperature that holds when e = £i is 



/3i = 



(i-q)rd^-i ■ 



(55) 



C. Approximate approach to the temperature 

The approach just taken in defining the initial condi- 
tions suggests a convenient way of obtaining an approxi- 
mate expression for the temperature. 

Consider "dilute" states j where cracks do not signif- 
icantly interact. In this case, the probability P™ that 
any one cell has broken when the maximum strain ten- 
sor is at e„i is again just the cumulative distribution 
P,n = [i^£,n ■■ SC : £^/{2Td^-^)]''+\ In this case, the 
probability of observing a non-interacting state j con- 
sisting of Afj cracks is pj — Pm^ [1 — PmY^~^^^ where 
we have forgone the analysis of the preceeding section 
demonstrating that the unbroken-cell probabilities are 
negligibly infiuenced by the strain perturbations from the 
Nj broken cells (at least for k > 0). We may write 



Pj = Po exp 



(56) 



where po = (1 ~ Pm)^ is the probability of the unbroken 
state. 

For such dilute states, the Hamiltonian of Eq. ( p6|) is 
written (with = e) 



E, 



1 



:Co:£„ 



(1-9) 



2 ^ 2 
so that our probability law predicts 
'f3{l-q) 



£„:(5C:£™7Vj (57) 



Pj =Po exp 



5C : SmNj 



(58) 
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Upon using 1/P™ = [2Vd°-^ / : SC : £mT^^ and 
equating Eqs. (|5^ and (|56|), the temperature is identified 



-2 In { [2rd^- V (£^£^ : SC : £„)]^-+^ - 1} 
(1 - q)em-SC:£m 



/?(£,„) - 

(59) 

This expression for has the expected behavior that 
/? = — oo when = 0, and that /3 is a negative and in- 
creasing function of up to the strain point Pm = 1/2 
where it smoothly goes to zero. For Pm > 1/2, /3 is a 
positive and increasing function of £„• Our probabihty 
law with (3 negative predicts the intact state to have the 
greatest probability, while when Pm > 1/2 and /3 is pos- 
itive, the most probable state jumps to every cell being 
broken. Although such a phase transition occurs in fiber 
bundles [Q, we demonstrate in Paper III using the ex- 
act differential equation for temperature, that the local- 
ization transition always occurs prior to this divergent- 
temperature transition. 



We emphasize that Eq. ( |59| ) is an approximation to 
the extent that due to the long-range nature of elastic 
interactions, one can never truly define a non- interacting 
state. We use it to obtain an order-of-magnitude idea of 
the temperature at a given strain. But it should always 
be considered preferable to obtain the temperature by 
integrating the exact Eq. ( p7| ) from the first-crack (or 
other exact) initial conditions. 



V. CONCLUSIONS 

The present theory of fracture in disordered solids 
works from the postulate that the probability pj of ob- 



serving a mesovolume in a given emergent-crack state j 
and at a given applied strain can be determined by maxi- 
mizing Shannon's measure of the emergent-crack disorder 
subject to constraints that come from the energy balance 
of brittle fracture. These constraints are what allow non- 
uniform probability distributions to occur. The vaHdity 
of this postulate can be demonstrated in simpler cases 
by integrating the probability distribution through 
history, but its general validity in the case of rocks with 
interacting cracks remains an open problem. Our ap- 
proach to answering this question is to use the statistical 
mechanics that follows from our maximal-disorder pos- 
tulate to make predictions about the physical properties 
of real systems and to compare such predictions to labo- 
ratory data. 
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